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Abstract 

Recent scanning tunnelling microscopy (STM) experiments on underdoped cuprates have dis- 
played modulations in the local electronic density of states which are centered on a Cu-O-Cu bond 
(Kohsaka et al, Science 315, 1380 (2007)). As a paradigm of the pinning of such bond-centered 
ordering in strongly correlated systems, we present the theory of valence bond solid (VBS) corre- 
lations near a single impurity in a square lattice antiferromagnet. The antiferromagnet is assumed 
to be in the vicinity of a quantum transition from a magnetically ordered Neel state to a spin-gap 
state with long-range VBS order. We identify two distinct classes of impurities: (i) local mod- 
ulation in the exchange constants, and (ii) a missing or additional spin, for which the impurity 
perturbation is represented by an uncompensated Berry phase. The "boundary" critical theory 
for these classes is developed: in the second class we find a "VBS pinwheel" around the impurity, 
accompanied by a suppression in the VBS susceptibility. Implications for numerical studies of 
quantum antiferromagnets and for STM experiments on the cuprates are noted. 



1 



I. INTRODUCTION 



A number of recent scanning tunnelling microscopy experiments have highlighted spatial 
modulations in the local density of states in the cuprate compounds, nucleated by external 
perturbations. In Ref. [l|, the spatial modulation was observed in the normal state above 
Tc, presumably nucleated by impurities. In Refs. the order was found in a halo 

around vortices, which were in turn pinned by impurities. Most recently, in Ref. js], similar 
charge-ordering patterns were found to be ubiquitous in the underdoped cuprates at low 
temperatures, and it was established that the charge ordering was "bond-centered", and 
had an anisotropic structure similar to a valence bond solid state^'^'^'^. 

In the light of these observations, it is of general interest to study the appearance of 
varieties of charge order (including "valence bond solid" (VBS) order-i^i^) near impurities 
in strongly correlated systems. For superfluid states, such a theory has been presented in 
earlier work^ii^, and compared quantitatively with some of the above experiments. It was 
argued that the charge order was hnked to quantum fluctuations of vortices/anti- vortices in 
the superfluid order. Consequently, the problem mapped onto the pinning of the vortices by 
impurities, and the quantum zero-point motion of vortices about the pinning site. In both 
zero and non-zero magnetic fields, enhanced charge order was found in the spatial region over 
which the vortex executed its zero-point motionii. This charge order was present even when 
the net vorticity was zero everywhere (as is the case in zero magnetic field): the vorticity 
cancelled between the vortex and anti- vortex fluctuations, but the charge order did not. 

This paper will present an extensive field-theoretic analysis of a paradigm of the problem 
of charge order near impurities in correlated systems. We will consider insulating 5* = 
1/2 antiferromagnets on the square lattice, across a quantum phase transition from the 
magnetically ordered Neel state, to a spin-gap valence bond solid (VBS) staie^'^^'^^^. By 
representing the 5 = 1/2 spins as hard-core bosons, our results can be reinterpreted as 
applying to the superfluid-insulator transition of bosons at half-filling on the square lattice: 
the Neel state of the antiferromagnet maps onto the superfluid state of the bosons, while the 
VBS state maps onto a Bose insulator with bond-centered charge order. The bond-centered 
charge correlations in the underdoped cuprates now appear to have two possible physical 
mechanisms ("disordered" antiferromagnet/superfluid), but it was argued in Ref. [l7| that 
they represent the same underlying physics. Our results here will go beyond the earlier 
workiii^ in two important respects: 

(i) We will describe the critical singularities in the impurity-induced VBS/charge order at 
the quantum critical point, and 

(a) We will consider a wider class of impurity perturbations. In the previous work— i^, an 
"impurity" was assumed to be a generic deformation of the underlying Hamiltonian which 
broke its space group symmetry. For the Neel- VBS transition, such an impurity is realized 
e.g. by the modulation in the magnitude of a particular exchange coupling - see Fig [H 
We briefly will discuss the critical singularities describing the enhancement of VBS order 
near such an impurity in Section II Al below; these results have a natural extension to the 
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FIG. 1: A modulated exchange impurity which is described in Section iLAl The dashed hne indicates 
a different value of the antiferromagnetic exchange constant. We expect VBS order to be enhanced 
near such an impurity, because the modulated exchange will lock in a preferred orientation and 
offset of the VBS state. 



models of charge order near the superfluid-insulator transition discussed above. However, 
the primary focus of the present paper is on a distinct class of impurities, in which the 
valence-bond structure of the non-magnetic ground state of the antiferromagnet is more 
strongly disrupted, and a "Berry phase" contribution of an unpaired spin is the crucial 
impurity- induced perturbatioiti^'^. Such impurities are realized by replacing the S = 1/2 
Cu spins in antiferromagnets by a non- magnetic Zn ion, or a S" = 1 Ni ion (see Fig. [2]). 
For the superfluid-insulator transition, such an impurity is a site from which particles are 
excluded, and so a local "phase-shift" is induced in the charge order of the insulating state 
(replacing a Cu atom by Zn or Ni is expected to have the desired "Cooper pair" exclusion 
effect^). Our main results will include a description of the suppression of VBS order near 
such "Berry phase" impurities: these results are summarized in Section II Bl below, and 
described in the body of the paper. A simple sketch of how such an impurity disrupts the 
VBS order is shown in Fig. [21 this figure builds upon the dual theory of spinous in the VBS 
state developed by Levin and Senthil^. The bulk of this paper will describe how quantum 
fluctuations of the type sketched in Fig. [2] lead to a modification of the scaling dimension of 
the VBS order in the vicinity of the vacancy. 

The remainder of the paper will be presented in the language of the Neel-VBS transition 
in quantum antiferromagnets. For this model, a field theoretic description of the vicinity of 
the quantum critical pointr^'^'^'^^ is provided by the CP^~^ theory at = 2: 



S = J (fxdr \{d^-iA^)z^\'^ + s\z^\^ + ^{\za\^Y + ^ie^^xduAxf 
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(a) (b) 

FIG. 2: A vacancy (the shaded circle) in a square lattice quantum antiferromagnet which is de- 
scribed in Section II Bl and the remainder of the paper. The thick lines represent singlet bonds 
between the spins. The local value of the VBS order is measured by the phase factors on the 
singlet bonds. Moving anti-clockwise from the right in both figures, we observe that VBS order 
cycles as 1 ^ i ^ — 1 — > — i in (a), and as i ^ — 1 ^ — i ^ 1 (b). Thus, both configurations are 
"vortices" in the VBS order, which we name "VBS pinwheels" (these VBS pinwheels are "dual" to 
the vortices in the superfluid/Neel order that are discussed in the beginning of the paper). Anti- 
pinwheels in the VBS order appear only around vacancies on the other sublattice; in other words, 
VBS pinwheels transform to VBS anti-pinwheels under translation by a single site — see Fig [3] later. 

Here /x, z/ A are spacetime indices, 2;^, a = 1 . . . is a complex scalar which is a SU(A^) fun- 
damental, and is a non-compact U(l) gauge field. The Neel order of the antiferromagnet 
is rf" = z'^T'^z, where T"' is a SU(A^) generator. The SU(A^) symmetry is spontaneously 
broken in the Neel phase, {n°') 7^ 0, which is realized for s < Sc, where Sc is the critical 
value of the tuning parameter, s, for the quantum phase transition. For s > Sc, the CP^^^ 
theory above describes a U(l) spin liquid state of the antiferromagnet, with gapped spinons 
Za and a gapless, U(l) photon. However, as has been argued at length elsewhere^^ii^, lat- 
tice effects not included in the continuum field theory (11.11) eventually render the U(l) spin 
liquid unstable to spinon confinement and fully gapped state with VBS order. The VBS 
order parameter, V, is an operator— which creates a Dirac monopole with total flux 27r 
in the U(l) gauge field A^. This paper will therefore be concerned with correlations of the 
monopole/VBS operator V under the field theory S after including the impurity perturba- 
tions described below. The bulk scaling dimension of the monopole operator at the s = Sc 
critical point will make frequent appearances in our analysis, and so we define this as 

= dim[V(x, r)] in the theory S without an impurity. (1-2) 
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The following subsections will now describe the two classes of impurity perturbations to 
the theory S shown in Figs. [1] and [2] respectively. 



A. Modulated exchange 

A modulation in the magnitude of an exchange constant in the underlying antiferromag- 
net (see Fig. [1]) breaks the lattice space group symmetry, but preserves the spin rotation 
symmetry. Also, the number of spins on each sublattice is preserved, so no "Berry phase" 
term is expected. Consequently, we need to consider all local perturbations to S which 
preserve the required symmetries. The simplest allowed possibility is a local shift in the 
position of the critical point. For an impurity at the spatial origin, a; = 0, this would lead 
to a term 



However, a simple computation^^ shows that 's is very likely an irrelevant perturbation at 
the bulk critical point. We have dim[s] = 1 — {D — l/u), where D = 3 is the spacetime 
dimension, and u is the correlation length exponent of S. Because it is almost certainly the 
case that u > 1/2, we conclude that s'is irrelevant. However, a more interesting perturbation 
is that considered in previous work— on the superfluid-insulator transition. In the present 
context, this perturbation follows from the fact that with broken space group symmetry, a 
linear coupling to the monopole operator is permitted. So we have the impurity action 



where his a complex- valued constant whose value depends upon the details of the modulated 
exchange near a; = 0. Now the renormalization group (RG) flow of h follows from Eq. (11. 2p 
to linear order 



The remainder of this subsection will analyze the correlations of the monopole/VBS operator 
V{x, t) in the theory S + iSimp. 

First, let us consider the likely possibility that < 1. In this case, h is a relevant 
perturbation, and higher order corrections to Eq. (II. 5p cannot be ignored. By analogy with 
results in the theory of boundary critical phenomenal, and in particular with the theory 
of the "extraordinary" transition^^'^'^, we conclude that a likely possibility is that the RG 
flow is to strong coupling, to a fixed point with \h\ = oo. In this, case some powerful 
statements on the correlations of V{x, r) can be immediately made. It is useful to express 
the correlations in the vicinity of the impurity by an operator product expansion (OPE). In 
general, this expansion will have the structure 






(1.4) 
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where l^mp is an operator localized on the impurity site, and Aj^p is the difference in scaling 
dimensions between V and Vimp- Specifically, Eq. fll.6p implies 



A^ = -ALp + dim[Mmp]. (1.7) 

Now at a \ h\ = oo fixed point, we expect that fluctuations of V near the impurity are strongly 
suppressed, and so it is a reasonable conclusion that Vimp is just the identity operator 

V^mp = 1. (1.8) 

Consequently, dim[V^inp] = 0, and we have our main result 

Ar^p = -A^- (1-9) 

The combination of Eq. (11.61) and (11.91) appears to be a promising route to measuring the 
scaling dimension of a monopole operator in numerical studies of quantum antiferromagnets. 

To complete our analysis of modulated exchange, we also address the case with A^ > 1. 
In this situation, by Eq. (II. 5p . the perturbation h is irrelevant, and so we may compute the 
consequences of h by perturbation theory. Computing correlations to first order in h we see 
that Eq. (II. 6p is now replaced by 

lim V(x, t) r^h 1x1"^'^''+^ (1.10) 

B. Missing spin 

Next we will consider the behavior of the monopole/ VBS operator V near the missing 



spin impurity illustrated in Fig. [2l As discussed in some detail in Ref. [iSj, the dominant 



consequent of such an impurity is an exactly marginal perturbation to S given by 

Simp=tQ [ dTA,{x = 0,T) (1.11) 



where Q is a "charge" characterizing the impurity. The value of Q does not flow under the 
RG, and so Q is a pure number which controls all universal characteristics of the impurity 
response. For an impurity of Fig. [2] with a single missing spin, Q = ±1. The remainder of 
this paper presents an analysis of the critical properties of the S + 5imp defined in Eqs. ( 11. ip 
and ffLTT]) . 

The magnetic correlations of the theory S + Si^p (and of a related theory^) have been 
computed in previous papers^^i^ which obtained the scaling dimensions of the Neel order 
parameter, n", and of the uniform magnetization density in the vicinity of the impurity. 
It was found that the impurity significantly enhanced the local magnetic susceptibilities. 
For the case of double-layer antiferromagnets, which have magnetic ordering transitions 
described by Landau-Ginzburg- Wilson theory, such impurity magnetic correlations have 
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FIG. 3: A VBS aniz-pinwheel in the presence of an impurity with a charge Q opposite in sign to 
that required by Eq. ()1.12p . This configuration has a higher energy cost than the VBS pinwheel 
configurations in Fig. [2j 



also been computed by similar methoda^'^'^, and found to be in excellent agreement with 
numerical studies^'^'^. 

This paper will describe the "charge-order" correlations of the theory S + 5imp by a 
computation of the OPE of the monopole/VBS operator V{x,t) as a; ^ 0. Our principal 
result is that the OPE is modified from the form in Eq. (11. 6p to 

lim V{x, t) ~ |f l^.'^p e-'^^ 14mp(r) (1.12) 

where 9 is the azimuthal angle of x. There are two important changes from Eq. (11. 6p . 
The first is that V^mp is no longer a trivial unit operator, but a fluctuating impurity degree 
of freedom with a non-trivial scaling dimension. The second is the presence of the e"**^^ 
factor, which indicates a Q-fold winding in the phase of the VBS order parameter around 
the impurity. The sketches in Fig. [2] give a simple physical interpretation of this winding 
in terms of the valence bond configurations of the underlying antiferromagnet: we will call 
this vortex-like winding in the VBS order a "VBS pinwheel." Also, as we discussed earlier—, 
the sign of Q is determined by the sublattice location of the missing spin. Thus, the result 
Eq. (I1.12P indicates that VBS pinwheels will occur preferentially around impurities on one 
sublattice, while VBS anti-pinwheels occur around impurities on the other sublattice. This 
same result is also obtained from the intuitive microscopic pictures in Fig. O Also, we show 
in Fig. [3] an illustration of an anti-pinwheel in the presence of an impurity on the disfavored 
sublattice: the same sublattice bond indicates that this configuration has a higher energy. 

Apart from establishing the form of Eq. (I1.12p . we will also describe computations of 
the exponent A^^p. There are general reasons for expecting that A^^p > 0, and this will 
be the case in the explicit result we obtain. This positive value of Aj^ characterizes the 
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suppression of VBS order near the impurity, and should be contrasted with the negative 
value in Eq. (11.91) for the impurity in Fig. [H 

Our analysis will begin in Section [III by a large analysis of the theory S + iSimp with 
full SU(A^) spin symmetry. We will establish Eq. (11.121) in this limit. We will also find that 
the N = oo limit (at fixed Q) of the exponent Aj^p vanishes, but we will not evaluate the 
subleading correction in the 1/N expansion here. 

The remainder of the paper will explore another approach to estimating A?J^p. This 
reliesi^ii^ on examining the "easy-plane" limit of the CP^~^ model, in which the global 
SU(A^) spin symmetry is reduced to f/(l)^^^. With this simplification to an abelian global 
symmetry, an explicit duality transformation of the theory becomes possible. In the dual 
theory, the monopole/VBS operator V has a local expression in terms of the dual fields, 
and so this facilitates the analysis of the impurity critical property. We will begin the dual 
analysis in Section IIIII by considering the simplest = 1 case—: this model describes the 
onset of VBS order in a = 1/2 quantum antiferromagnet in the presence of a staggered 
magnetic field^, and is the simplest setting in which several technical issues can be described. 
We then extend the analysis to general N in Section HVl The exponent will be estimated 
in these sections by a self-consistent theory of Gaussian fluctuations about a mean-field state; 
in the physically interesting case of = 2 and Q = 1, which describes both the easy plane 
antiferromagnet and the boson superfluid/insulator transition, we obtain the estimate 

ALp^O.57, N = 2,Q=1 (1.13) 

Our analysis of the easy plane theory in Section [IV] also exhibits certain features which we 
do not expect to be shared by the case with global SU(A^) symmetry: for Q/N = 1/2, we 
find VBS-vortex solutions in which the e"*'^^ factor in Eq. (11.121) is replaced by e~*^^ with 
the integer —Q < i < Q. In the self-consistent theory we present here, all the values of i 
are degenerate, but we expect these degeneracies are partially lifted in the full easy-plane 
theory. These issues are discussed further in Section IIVI and in a forthcoming paper. 



II. 1/N EXPANSION OF THE CP^"^ 
MONOPOLES 



THEORY IN THE PRESENCE OF 



The insertion of one monopole into the partition function of CP model in the dis- 

The 1/A^ expansion proceeds by 
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ordered phase has been originally considered in Ref. 
replacing the quartic self-interactions in Eq. (II. ip by a fixed-length constraint on the spinous; 
so we consider the action 



S 



d xdr 



\{d^-iA^)za\'^ + iX \z, 



(2.1) 



where A is a fluctuating Lagrange multiplier field. The procedure for generating the 1/A^ 
expansion is now simple. One first integrates over the z fields obtaining an effective action 
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for Afj^ and A. However, instead of expanding this effective action around the trivial classical 
vacuum = 0, one expands around the monopole (instanton) solution, A^, with 



= 2nq^ — (2.2) 

where = e^u\dijA\, q is the monopole charge and xq is the monopole position. In practice, 
integrating out the z fields in the background of spatially varying monopole fields is quite 
complicated (even more so due to the appearance of UV and IR divergences), so that only 
the leading term in the expansion has been computed in the past (that is fluctuations 
of Afj^ about the monopole solution have not been taken into account). At this order, one 
finds. 



{V'{x)) ~ [j) (2.3) 

where V(x) is the monopole operator of charge q, m is the mass gap of the theory, A is the 
ultraviolet cut-off and pq is a collection of universal numbers (depending only on the charge 



of the monopole) which have been computed in Ref. |2l|. Thus, the dimension of operator 
Vi{x), dim[\/'?] = 2Npq. 

If finding the expectation value of a monopole operator (and its scaling dimension) was 
very complicated, finding correlators of V{x) with Wilson loops at = oo turns out to be 
exceedingly simple. Indeed, we notice that at leading order in 1/A^ it is sufficient to simply 
replace in the Wilson loop by its monopole value. 



(V«(x)exp [-iQ j^A^dx^)) 



exp (^-iQ A^dx^ = exp (^-iQ F'/S^ (2.4) 



provided that we take the charge Q to be 0{1) in A^ (otherwise, if Q ~ 0{N) the Wilson line 
will change the background monopole field and the problem becomes intractable). Here C is 
some closed contour and S is any surface such that dS = C. Thus, all we have to do is find 
the flux of our monopole through the Wilson loop that we are considering. Fluctuations of 
Afj^ about the monopole field (12. 2p will contribute at 0{1/N) to the correlator (12.41) . Likewise, 
if we denote the Wilson loop operator by W{C), then in the absence of the monopole field 
{W{C)) ~ 1 + 0(1/A^) (saturated by fluctuations of A^ around the trivial vacuum), so 



{V'^{x)W{C)) 
{W{C)) 



{V'^)exp(^-tQ j^F^dS,^ (2.5) 



and at leading order in 1/A^ the external charge only changes the phase of the expectation 
value of monopole operator but not its magnitude. 

In principle we are interested in finding the correlator of the monopole operator (that 
we place at a point x = {rcos6,rsm6,0)) and a straight, temporal Wilson line of charge 
Q (which we place at the origin). However, to regularize possible IR divergences let's also 
place a charge —Q on the positive x axis far away from the origin. As usual, we may connect 
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the two oppositely directed Wilson lines in the far past and far future. Then, according to 
fl2.4p we have to compute the magnetic flux due to the monopole field (12.21) through the 
y = 0,x > half-plane, 

F-dS = -I r dr r dx rsin^_^ _ ^^.g) 

2 J-oo Jo ((x — rcos6')2 + r^sin ^-|-t^)2 

oo 



^^7—7777^^^^:7—27 (2-7) 



r sin 9 

(x — r cos Oy + sin^ 6 
-qin - 9) (2.8) 







We see that the flux through the Wilson loop changes by 2Txq as the monopole crosses the 
surface of the loop. However, the expectation value, 

(^^(a^))imp = ^^''l^^c))^^ = (^^)e^««(-^) (2.9) 

remains single valued, as by Dirac's condition Q is an integer (in what follows, we shall also 
often discuss Wilson loops with non-integer charge Q, which in the presence of monopole 
operators are defined by specifying a surface S, W{S) = e~^^ ^s^t^'^^i^. The correlation 
functions then explicitely depend on the choice of the surface, as can be seen from (12.91) ). 
Thus, we see that the phase of the monopole operator winds by —2ttQ as we move it in a 
full circle around the Wilson line, i.e. an external charge creates a vortex of the monopole 
field, consistent with the OPE in Eq. I1.12[ We expect that once we go beyond the leading 
order in A^, this vortex will also get a nontrivial spatial profile, 

(^'(x))imp = (y^)f{m\x\)e-''^^'e'^ (2.10) 

Here /(r) is the vortex profile function and e*^ is some overall phase (discussed below). 
We expect that far away from the external charge, the monopole field tends to its vacuum 
expectation value so that f{oo) = 1. Moreover, by continuity we expect the monopole field 
to vanish at the origin, /(O) = 0. To the order to which we were working, /(r) = 1, which 
imphes that the impurity exponent Aj^p ~ 0{1/N). 

Notice that the result ( 12. 9p is sensitive to the angular position of the distant charge 
relative to the one at the origin (we introduced the variable 6 as the angle between the plane 
of the Wilson loop and the monopole operator). This is not unexpected: the monopole 
field is the order parameter for the flux symmetry, which is spontaneously broken in the 
disordered phase. As we rotate the distant charge, the overall phase e*^ of the expectation 
value of the monopole operator changes - that is we explore different states in our vacuum 
manifold. 

If we were instead considering a correlation function of a string of monopole operators 
Yii V"'^'(xj) such that the overall combination is invariant under the flux symmetry (that is 
— 0) ^6 expect the dependence on the angular position of the distant charge to drop 
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out. We can check this in the hmit m\xi — Xj\ ^ 1, m\xi\ S> 1, assuming clustering, 

(n^''(^^))imp - n<^''(^^))'-p ^ ii^vf)e-^^^^'^ (2.11) 

i i i 

which is invariant under 6i — * 0i + X- Alternatively, in the same limit of far separated 
monopoles and at = oo, the classical magnetic field will just be a linear superposition 
of magnetic fields due to each monopole. Thus, the flux $ through the Wilson loop will 
be given by, $ = — ^^(vr — 6i) = qi9i and using the equivalent of (12. Sp for a string of 
monopole operators, we arive at the same expression (12. lip . 

We expect the general form (12.101) to be preserved at any finite order in 1/A^. Nevertheless, 
in the fiux-broken phase of the theory, there are also non-perturbative effects that should 
be taken into consideration. Indeed, the ?7(1)$ vortex nucleated by the external charge is 
global, and thus, will have a logarithmically divergent energy. Put into a more conventional 
language, the external charge creates a Coulomb potential, which is logarithmic in two 
dimensions, V{r) ^ —^^\og{mr) for mr ^ 1. The effective coupling constant can be 
calculated in the 1/A^ expansion to be ~ jj^- Thus, it will be energetically favourable 
for the external charge to bind a dynamical spinon (we concentrate on the case Q = 1 here 
for simplicity). This process can be analyzed by means of a non-relativistic Schrodinger 
equation^. One finds a bound state of size ~ Nim~^. We expect that for r ^ r-f, the 
external charge will be screened by the dynamical spinon. On the other hand for r 
this logarithmic confinement should generally have little effect on the physics. However, 
there is one notable exception: the expectation value of the monopole operator (12. 9p 
will be drastically altered on all distance scales by the screening. Indeed, if we assume that 
screening takes place, {V^l^x)) has to tend to its vacuum expectation value for ^ r^, 
and should experience no phase winding. We don't expect the winding number to change 
abruptly as we decrease \x\, so we won't see a phase winding of {V^i^x)) on short distances 
\x\ <^ Tf, as well. 

A toy model for the disappearance of winding when screening effects are taken into 
account can be constructed as follows. We can use the charge —Q that we previously put 
far away from the origin to represent the dynamical spinon that gets bound to the external 
charge. We first freeze the location of this spinon at some position x' away from the origin 
and compute the resulting expectation value of V^i^x) using eq. ( 12. Sp . We then average the 
resulting (V^'^(x)) over the spinon positions x' with the probability distribution |'?/'(x') p, where 
iplx) is the spinon wave-function. Since this wave-function will be azimuthally symmetric, 
one immediately learns that upon averaging over the angular position of the spinon, (V^lx)) 
looses its finite winding number and will, in fact, carry a constant phase for all x. This 
same averaging will also lead to an additional supression {V^i^x)) ~ \x\ as x (recall 
that at A^ = oo there was no supression of the vortex profile for x ^ before screening 
effects were taken into account). The origin of this supression is easy to see - for an external 
charge located infinitely far away, the averaging over the azimuthal position of the charge is 
identical to averaging of the phase x sq. (I2.10p producing a zero result for (V^). 
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Do the above findings invalidate the OPE fll.l2p ? The answer is no. The above discussion 
simply implies that (V^mp) = and, thus, the expectation value {V{x)) ioi x ^ is controlled 
by higher order terms in the OPE (namely, by the impurity operator with angular momentum 
zero). However, higher correlation functions of V operator, e.g. the VBS susceptibility 
(V^(x)K^(x')), are still controlled by the OPE (11.12p . Such correlators are invariant under 
the symmetry, so as we argued above, their short distance properties are not sensitive 

to the location of the distant charge, and hence, to screening physics. 

III. EASY PLANE MODEL AT = 1 

This section, and the next, will examine a simplified version of the theory S + iSimp in 
which the non-abelian global SU(A^) symmetry is reduced to an abelian U(l)^~^ symmetry. 
This enables us to use the tools of abelian particle-vortex duality^i^ to obtain a theory 
expressed in terms of fields which are locally related to the monopole/VBS operator V. The 
present section will consider the simplest case^ with = 1. This model describes the 
onset of VBS order in a S* = 1/2 quantum antiferromagnet in the presence of a staggered 
magnetic field^, and is useful in resolving a number of key technical questions in their 
simplest setting. For = 1, the theory S does not have any global continuous symmetry, 
and becomes equivalent to scalar electrodynamics. With the results for the = 1 theory 
obtained in the present section, we will be able to rapidly analyze the general A^ case in the 
next section. 

A. Duality and Wilson loops 

It is well known that in three space-time dimensions, near its critical point, non-compact 
A^ = 1 scalar electrodynamics is dual to a theory of a complex (pseudo) scalar field with a 
global U{1) symmetry2^i^. The Lagrangians of these two theories are as follows, 

Lqed = ^F'^ + {{d, - tA,)z\' + m'\z\' + ^Izl"" (3.1) 
LxY = \d,V\' + m'\V\' + ^\V\' (3.2) 

Here z and V are complex one component fields. The duality is understood as being true 
for the range of parameters where Lqed has a second order phase transition (which at weak 
coupling is believed to occur for g/e'^ sufficiently large). One way to understand the duality 
is by noting that the phase transition in scalar QED is driven by spontaneous breaking of 
flux symmetry f/(l)$, which is precisely the global symmetry of Lxy- The order parameter 
for the flux symmery is the monopole operator V{x) - that is the dynamical field of Lxy^- 
As we know, to each continous symmetry there corresponds a conserved current. In the case 
of flux symmetry of QED, this pseudo-vector current is just the magnetic field F^, which is 
trivially conserved in the absence of monopoles, d^Fn = 0. Let's introduce an external field 
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iJ^j that would couple to this current, 

^Lqed = iH^Ff, (3.3) 

Suppose we are calculating some correlation function with insertion of a string of monopole 
operators of charge at points Xi. The gauge field in the path integral is then subject 
to the condition, dfj,F^ = ^j27rgj5(x — Xi). Then under the transformation, 

H^^H, + d,a (3.4) 

Sqed Sqed + i j dx d^aF^ = Sqed -i J dxa d^F^ = Sqed - Svrz ^ qia{xi) (3.5) 

Hence, by introducing the field we can enlarge the global t/(l)$ symmetry to a fictitious 
local symmetry, provided that the monopole operators transform as, 

Vi{x) ^ e'"'^"("V''(a;) (3.6) 

The dual Lagrangian Lxy has to posses this local symmetry. Hence, to introduce the field 
Hn into the dual Lagrangian we simply have to covariantize the derivative of the dynamical 
monopole field V, 

d,V ^ D^V = {d^ - 2mH^)V (3.7) 

in eq. (13.21) . Other "gauge invariant" operators can also be added to Lxy, e.g. H^i,; however, 
their contribution will, generally, either cancel out in correlation functions or be less singular 
near the critical point. 

Thus, the dual Lagrangian in the presence of a background source field if^ is given by, 

Lxy =\{d^- 2mH^)V\^ + m^\V\^ + ^\V\^ (3.8) 



The covariantization procedure (13.70 was explicitly written down in Ref. l38|. Similar 
arguments for the case of a constant imaginary iJ^, which physically represents an external 
magnetic field in the QED language and translates into a chemical potential for the flux 
symmetry in the XY language, have been given in Ref. |39|. In a companion paper— we shall 



also give an argument based on an exact duality transformation on the lattice, which will 
support (13.81) . 

Having learned how to incorporate the source field if^ into the dual Lagrangian, it is now 
trivial to dualize Wilson loops. Indeed, insertion of a Wilson loop W{C) into a correlation 
function is equivalent to adding into the Lagrangian the source term 

5L = iQ J dXf^Af, = J dS^F^ = i J dxH^^F^ (3.9) 
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where 

H^ix) = Q [ dS^5{x-y) (3.10) 

Jy&S 

That is is a field that fives on the surface of the Wilson loop and is directed perpendicular 
to this surface. 

Another benefit of introducing the source field is that by differentiating with respect 
to it we can compute correlation functions of the magnetic field F^. For instance, 

{-^F,{x))h = = -2^^{[V^D,V - (D.Vyv) {x))h (3.11) 

That is the topological fiux current of QED gets mapped into the Noether's current 
associated with the global U{1) symmetry of the dual model. Differentiating once again, 

(F,(x)F.(,)),,_ ^ -jj^^^^ (3.12) 



= {271 f (^(0 D ^V{x)V^ D .V{y))H,conn + 26^Jix - y) {VW{x))h^ 

(3.13) 

The first term in fl3.13p is the expected correlator of two U{l)q, currents, while the second 
term is a tadpole that ensures the overall transversality of the correlation function. 

Having discussed the duality at length, we now return to our original problem: what 
is the infiuence of the exernal charge (Wilson line) on various physical observables. The 
observable of most interest to us is the monopole operator V{x). However, this observable 
is physical only for integer- valued charge Q of the Wilson line (Dirac's condition). Indeed, 
recall that in the dual language the field H depends on a choice of surface S of the Wilson 
loop. If we pick a different surface S' then the field undergoes a gauge transformation 
H'^ = + d^a with a{x) = —Q l^^gv where V is the volume bounded by the two 
surfaces S and S'. Hence, 

(l^(x)...)H' = e2™(^)(r(x)...)^ (3.14) 

where ellipses denote some other operators. Thus, the operator V{x) is invariant under 
changing the surface of the Wilson loop if and only if Q is an integer. However, if the charge 
Q is a rational number, Q = p/q where p and q are integers then the fiux 27rg monopole 
operator V'^{x) ~ {V{x)y is physical. Moreover, a theory with arbitrary irrational Q is still 
sensible provided that we confine our attention to correlation functions of operators which 
are invariant under the fictitious f/(l)$ local symmetry, e.g. the magnetic field operator 
—iFf^ = In fact, if we are dealing with such gauge invariant operators we 

don't necessarily have to use the precise form of H given by fl3.10p : defining 7^ to be a field 
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living on the perimeter of the Wilson loop and directed along it, 



we see that 



lpi{.x) = Qj dy^6{x-y) (3.15) 

lyGC 



^^,uX^uHx = 7m (3-16) 



Then, by performing a suitable gauge transformation on if^ and V we can choose if^ to be 
any field with curl given by 7^. Thus, we see that the duality maps a Wilson loop of charge 
Q in the QED language to an external magnetic flux tube of flux 2nQ in the XY language. 



This correspondence has been noted in Ref. |4l|, but the consequences of this correspondence 



for the critical properties of Wilson loops were not discussed. 

Now we can address the problem that we originally posed in a dual language. Let's place 
an charge external charge Q at the spatial origin. For now we don't insist that this charge 
be an integer. The dual source field if^ must, therefore, satisfy 

V xH = Q6^{x)f (3.17) 

Thus, we basically have to solve an Aharonov-Bohm problem with flux 27tQ. One choice for 
the source field if^ is 

H^{x)=Q6^,20{x)6{y) (3.18) 

This is the so-called string gauge, which corresponds to (13.101) . with the surface of the Wilson 
loop being the plane y = 0, x > 0. As is well known, the string gauge is equivalent to = 
and the boundary condition, 

V{e = 27t) = e-^^'^V{e = 0) (3.19) 

where 6 is the azimuthal angle. Thus, we have to solve the theory (13. 2p with the twisted 
boundary condition (I3.19p . We observe that the physics is, therefore, a periodic function of 
Q. For integer Q the boundary condition (I3.19P is trivial - there is no twist. So our argument 
indicates that integral external charges do not affect correlation functions on distances of 
order of the correlation length of the theory: screening of integral charges takes place on 
distance scales of oder of microscopic UV cutoff. This surprising fact is discussed in more 
detail in a companion paper—. 

The behaviour at non-integer Q is less unexpected. One physical question that we may 
ask is what is the magnetic (electric) field induced by the charge Q (we define the electric 
field Ei = = —eijFj where latin letters j, k run over spatial indices). Although this is 
a departure from our original goal, we will see that a lot of the results that we will obtain 
along the way will be useful when we return to discuss correlators of monopole field for the 
planar theory with N fields. Another question that we will adress for non-integer, rational, 
values oi Q = p/qis the behaviour of higher flux monopole operators V{xy. 
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B. Perturbative expansion of the dual theory for Q — > 

The magnetic field — iF^ is a conserved current and receives no renormalizations and, 
thus, has conformal dimension 2. Therefore, at the critical point we expect, 

{-iE)=C{Q)^f (3.20) 

The electric field is imaginary as we are working in Euclidean space. The coeficient C{Q) 
is a universal number that is a periodic function of charge Q. We shall be interested in 
determining this function. 

For Q — we can perform a perturbative expansion in ~ C>{Q). 

(3.21) 

As we have learned, the correlation function of magnetic field dualizes to, 

K,A^ -y) = {F,{x)FM) = (2vr)' {iy%V {x)V%V {y)) + 25,J{x - yW^)) (3.22) 
By transversahty, 

K,.(p) = ir(p)(V-^) (3.23) 
By RG K{p) should have the form, 

K{p) = Mg{p/M) (3.24) 

where M is some physical scale in the theory (e.g. in the f/(l)<i> disordered phase, the mass 
of the monopole field V). At the critical point, 

K{p) = A\p\ (3.25) 

where A is some universal number. On the XY side of the theory, this universal number 
has been computed before using both e expansion^ and large M expansion--. The large M 
expansion is obtained by replacing the complex scalar V in the action for the XY theory 
(13. 2p by an M component complex field. In the large M expansion the coefficient A is found 
to be at next to leading order in M, 

, / 1 32 \ M=i 

while in the e expansion one obtains A ^ 2.0 at 0{e^). Monte-Carlo simulations on the XY 
model^ indicate A ^ 1.8.— The coefficient A can also be computed by performing a large 
N expansion in the original QED, whereby the field z is promoted to have components. 
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At leading order one obtains A = 16/N 16 (as usual, direct large N expansion in QED 
produces results, which are numerically notoriously inaccurate for ~ 1). 

For completeness, we also discuss the behaviour of K{p) at small momenta on both sides 
of the critical point. In the phase where the f/(l)$ symmetry is spontaneously broken the 
spectrum of the theory should contain a goldstone, which can be created out of the vacuum 
by the U(l)q, current, 

lim(p|F„(a;)|0) = 27rlim(p| - iV^K^i^M = '^rrifp.e'P'' (3.27) 

p— >o p-^O 

where in three dimensions defines a physical energy scale. Note that equation (13.271) is 
written in Minkowski space. We see that the goldstone is nothing but the photon of the 
original QED. Then K^i,{p) should have a pole at = and using spectral decomposition, 

hmK{p) = {27rff (3.28) 

On the other hand, in the the phase where the U{l)is, symmetry is unbroken (that is in the 
"superconducting" phase of QED) the V field is massive and all the excitatations have a 
gap. Therefore, K^j,^{p) cannot have a pole at = and 

li^Kip) ~ |1 (3.29) 

Having discussed the expected form of K^^ in different phases we can go back to eq. (I3.2ip 
for electric field induced by the charge Q. Introducing the kernel 'D{p) = K{p)/p'^, and using 
eq. fl3T7D . 

i-tF^ix)) = - y dyK^^ix - y)H,{y) = -Q j dr'e^^^d^Vix, r') (3.30) 

Hence, 

{-iE{x)) = Qh{\x\)f (3.31) 

where 

/^(|x|) = -^y" dT'V{x,T') (3.32) 
Substituting the expression (13.251) for K{p) at the critical point we obtain, 

(-zE(f)) = Q^A_f (3.33) 

Hence we identify, 

C{Q) ^ QA/{27r), (3.34) 
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Similarly, in the ?7(l)<i> ordered phase, 

i-iEix)) = Q^^r (3.35) 

\x\ 

So in this phase, as expected, the external electric charge produces the usual Coulomb-like 
electric field, E = "^^^ , as appropriate to two spatial dimensions with the identification 

eeff = 271 f. 

C. Peculiarities of the free theory 

So far we have only discussed the leading term in C{Q) for Q 0. In principle, we could 
continue the expansion in Q to higher orders: then the problem reduces to finding correlators 
of current operators —iF^ = iV'^ 9^ V . These correlators can be found by performing either 
e or 1/M expansion of the XY model. In either case, going beyond the leading order in Q 
is not simple. So, instead, we choose to return to the formulation of the problem involving 
the twisted boundary condition (13.191) . In the next section we will use this formulation to 
compute C{Q) for all Q (albeit numerically) at M = oo. However, before we do so, we will 
solve a slightly simpler problem: namely we find the form of C{Q) at the gaussian fixed 
point ^ = 0, fh? = of the Lagrangian (13. 2p . The reason for studying the free theory is that 
the calculations in it are, technically, very similar to those in the strongly coupled M = oo 
theory addressed in the next section (even though the physical results are quite different). 

In the free theory, C{Q) can be determined exactly, and, surprisingly, turns out to be 
a non-analytic function of Q at Q = 0. We have not been able to see any hints of this 
non-analyticity from the perturbative expansion of the free theory in Q (perhaps because 
we could go perturbatively only to linear order in Q, whereas the non-analyticity of C{Q) 
starts only at order |Qp). On the other hand, once we go in the next section to the strongly 
interacting fixed point (obtained in the M = oo limit), the theory cures itself of all IR 
divergences and C{Q) becomes analytic in Q. 

So, let's compute, 

{-iF^{x)) = {-2tiiV%V{x)) = -27^1 lim(9^ - dl){V{x)V\y)) (3.36) 

in the free theory, L = subject to boundary condition (I3.19p . As eq. (I3.36P shows, to 

find the f/(l)$ current it is sufficient to determine the propagator, D{x — y) = {V {x)V'' (y)) . 
The propagator will also determine the correlation function of operators {V{x)y for rational 

Q ~ p/ 

{iVix)r{V\y)y) = q\Dix-yy (3.37) 
We note that our problem is invariant under translations along the temporal direction, so, 

D{x, x',t-t')= f ^D2{x, x', cu2)e'"("-"') (3.38) 

J 27^ 
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where D2{x, x', u"^) denotes the two-dimensional propagator with mass = u'^ and twisted 
b.c. (13.191) . We use spectral decomposition to find D2, 



^ gi(y poo 1 

D,{x,x',m') = J2ir dE—^—=ct>i,E{mA^) (3-39) 
^ 271 Jo + h 

where we sum over states with fixed azimuthal angular momentum / = n — Q, n G Z. 
Note that the angular momenta are not integral due to the twisted b.c. (13.191) . The radial 
eigenfunctions (pi, sir) satisfy, 



l§-,ir§-^) + lr.]Mr) = EMr) (3.40) 



and are normalized as, 



POO 

/ c?rr0*^(r)0z,s,(r) = 5{E - E') (3.41) 
Jo 



The solution to ODE fl3:i0|) is. 



hEir) = ^JmiVEr) (3.42) 



where Jn{u) is the n-th order Bessel function. Hence, 

fl(r,r',9-9',r-r') = $:e"«'-»') / ^e-"-'' |^J|„(.r) J„|(.r') (3.43) 

where we made the substitution u = \/E. Integrating over a;, 

1 poo I I 

Diry,e,T) = ^J2'''l dvJ\,\{^v)J\^{v)eM-^-^v) (3.44) 

Now we can ask, what is the behaviour of the propagator D{r,r' ,6,t) for r — >■ 0, i.e. for 
r < r'. Recalling, J\i\{r) ^ 2i' l r(K|+i) ^'''' 

dvJ\i\i^v)J\i\i^)^M-yv) ^ (3.45) 
with 

Thus, for r ^ the contribution of states with angular momentum / to the propagator scales 
as rl'l. So, the largest contribution comes from smallest |/| = \n — Q\. For — | < Q < ^ 
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smallest |/| is given by setting n = 0, I = —Q. Hence, for |Q| < 1/2, and r/r' ^ 1, 

1 / r \ IQI T 

(3-47) 

For values of |Q| > 1/2 we simply periodize the eq. (13.471) . since all physics in XY model is 
periodic in Q with period 1 (see discussion in previous section). From here on, we therefore 
confine our attention to \Q\ < 1/2. 

Thus, if we were to perform the OPE in Eq. (11.121) in the XY model 

V{x,t) ~ \x\^^^^e-'^^Vi^p{T) for \x\ (3.48) 

we would obtain for |(5| < | in the free XY model, 

ALp=IQI- (3.49) 

We immediately see that the free theory is non-analytic in Q at Q = 0. By periodizing in 
Q, we also see that Aj^p is non-analytic at Q = ±1/2. However, this later non-analyticity 
appears only after we take r — > limit of the propagator, while we expect the non-analyticity 
at Q = to persist in the propagator for arbitrary r, r'. 

In fact, Q = I is a very special point. At this point the n = 0, I = —Q and n = 1, 
I = 1 — Q, i.e. I = ±1/2 terms in the sum (13.441) become equally important for r/r' —>■ 0. 
Thus, for Q — > 1/2 it makes sense to keep both terms in the assymptotic expansion of the 
propagator, 

« 5^ ((^)'e-«%(^) + (^)"''e-»-""B«-(^)) (3.50) 

and we may hypothesize the impurity OPE, for Q — > 1/2, 

V{x,t) ~ CQ|f|^Se-«^rQ(r) +CQ_i|f|^S-ie-(«-i)^ VQ_i(r), for |f| (3.51) 

where Vq and Vq_i are two impurity operators, with impurity anomalous dimensions Aq 
and Ag_p In the free theory, = Q and A^_^ = 1-Q. Hence, for Q < 1/2, < Ag_^ 
and the operator Vq is the most relevant as |x| — 0, while the operator Vq_i provides a 
subleading correction. For Q > 1/2 the roles of these two operators are reversed. Finally, 
for (5 = 1/2 the two operators have degenerate anomalous dimensions, A^^^ = A^^^^ and, 

V{x, t) ~ ci/2|f|^V2e-^^/2 Vy2{T) + c_i/2|f|^-V2e''/2 V^-i/2(r), for |f | ^ (3.52) 

Physically, the Q = 1/2 point is special because the CP symmetry is effectively restored 
at it^. Indeed, under CP, Q —Q. However, as already discussed, the universal physics 
is periodic in Q, so the points Q = ±1/2 are identified. Thus, the two impurity operators, 
V±i/2 are just CP conjugates of each other and must have the same impurity anomalous 
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dimensions. Hence, although our original analysis was performed for the case of the free 
theory, we expect the conclusions to remain valid in the strongly interacting theory. 

We remind the reader that even though the operator V{x) is mathematically well defined 
by specifying the surface S of the Wilson loop for arbitrary Q, it is not physical for non- 
integral Q. Indeed, a physical operator cannot obey twisted boundary conditions. However, 
for rational Q = p/q, the flux 27rg monopole operator V^lx) ~ {V{x)y is well-defined on 
both sides of the duality. Using fl3.37p and fl3.47p . we obtain the OPE, 

r) ~ \x\^L,('})e-"iQ<^ l^^p(r) for |x| ^ (3.53) 

with 

ALp(9) = <l\Q\ (3.54) 

in the free XY theory for \Q\ < 1/2. Since qQ = p is an integer, the OPE fl3.53p is invariant 
under 6 6 + 2tt, making the operator \^^(x) single- valued, as required. 

Having discussed the impurity OPEs, let us return to the calculation of electric field. 
Since we know that the electric field will be radial, we only need the 6 component of the 
magnetic field, 

/-iFe) = -2m- lim (dg - dg') Dir = r',e- 9', t = r') = -Am- lim deDir = r', 9,t = r') 

r 0^9' r 9^0 

(3.55) 

For this purpose, we don't need the propagator with r/r' <^ 1, but rather with r — > r', 
r r'. We denote, D{r, 9) = D{r = r', 9,t = r'). Unfortunately, if we plug r = r', t = t' 
into the expression for propagator (13.440 . the integral over v diverges. We expect that if 
we instead first keep r — r', t — t' finite, perform the integration over v, sum over angular 
momenta I and only then take r = r', r = r', the divergence disappears. There are also 
other ways to regularize the propagator: e.g. make the integral over u in (I3.38P run over 
D — 2 dimensions. This would correspond to the XY model in D dimensions coupled to an 
external fiux-tube (the fiux-tube is a defect in 2 dimensions, so its world-volume is — 2 
dimensional). One then takes the limit — > 3 at the end of the calculation. We have 
successfully used this method to compute the electric field (see Appendix . The result for 
the coefficient C{Q) of eq. (Km is, 

C(Q) = ^(l-2|Q|)2tan(7rQ), |Q| < 1 (3.56) 

o 

Thus, we see that the function C{Q) is non-analytic aX Q = 0. This analyticity occurs at 
non- leading order in Q, 

C{Q)^^Q{1-A\Q\), Q^O (3.57) 

The leading order term, C (Q) ^ | Q is the one which would have been predicted by ex- 
panding the free theory perturbatively in Q. 

One can also derive the result fl3.56p in a different way, which can be more easily gen- 
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eralized from the free theory to the 1/M expansion in a strongly interacting theory. This 
calculation is based on the integral representation of the propagator of the twisted theory 



derived in Ref. |47|. We repeat the calculations of Ref. |47| in Appendix [B] as in the next 



section we will need to generalize them for application in 1/M expansion. The result is, 
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D{r,9) = —I dutanh{7Tu)U^{9) (3.58) 



with, 

^ e-^-Q-g'^W sinh(zy|g|) + sinh(z/(27r - \e\)) 
''^ ' cosh(27rz/) - cos(27rg) ^ ' ' 

from which one recovers eq. (I3.56P by using eq. (I3.55p . see Appendix [Bl 



D. 1/M expansion of the dual theory 

We now progress from the free XY model to the 1/M expansion of the strongly interacting 
theory. We take the Lagrangian to be, 

L=\d^V\^^%\{\Vf--^ (3.60) 

Here is an M component complex scalar and A is a Lagrange multiplier, which enforces 
the local constraint, 

|V^r = - (3.61) 
9 

This hard constraint replaces the self-interaction of the V field. In the presence of an external 
charge in the direct theory, we take V to satisfy the twisted boundary conditions (13.191) . In 
principle, we would like to solve the theory (13.601) in the limit M 1. However, practically 
we will only be able to perform computations at M = oo. 

We will be interested in the properties of the theory (13.191) at its critical point g = Qc- 
As is well known from standard 1/M expansion techniques, at M = oo the critical coupling 
is given by, 

^iV^V) = D{x = x') (3.62) 



Mg, M 

where D is the usual massless 3D propagator, 

D{x,x') = — --^ (3.63) 

Of course, the propagator with x = x' m. (I3.62p is UV singular and has to be regularized. 
Since we will perform calculations of propagator in position space, it is convenient for us to 
use point-splitting regularization. 

In the absence of the twisted boundary condition (13.191) and at the critical point, we 
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perform the expansion around (iX) = (so that the effective mass for the V particles 
vanishes). However, once Q is finite, A = is no longer sufficient to make the constraint 
f l3.6ip satisfied. Instead, the Lagrange muliplier aquires a spatial dependence 

(^A(f,r)) = ^ (3.64) 

Here a is a universal function of the charge Q. The dependence on x is determined from the 
canonical dimension of A (A aquires a non-trivial anomalous dimension only at order 1/M). 
Thus, at finite Q, the propagator of V field satisfies, 

i-d^ + j^)Dix, x', Q) = 5{x - x') (3.65) 

and a{Q) should be determined self-consistently from the equation. 

Combining eqs. fl3.62p . (13.661) . 

lim (L)(x, X , Q) - D{x, x' , Q = 0)) = (3.67) 

Thus, the problem is reduced to finding the propagator D{x,x',Q). Just as in the free 
case, we use spectal decomposition (13.391) . and the radial functions (pi,E{r) now satisfy, 

where again due to the twisted boundary conditions I = n — Q, n E Z. The solution to 
dMHD is, 

Mr) = ^J^(v^r) (3.69) 

Comparing the result above to free theory (I3.42p . we see that the only difference is in 
the replacement of the indices of Bessel functions |/| —>■ ^/P + a. Going from 2D to 3D 
propagator as in the free case (I3.44p . 

1 /"°° r Irl 

D{ry,e,r) = ^E^'' / dvJ^^{-v)J^^{v)eM'^-^v) (3.70) 

Finally, expanding the propagator ( ]3.70p for r ^ r', we obtain the equivalent of ( ]3.47p . 



1 / r \ \/Q'^+a.(Q) T 



23 



Thus, we recover the OPE (13.481) . but the impurity exponent now becomes some nontrivial 
function of Q, 

ALp = vW^, IQI<l/2 (3.72) 

We note that, similar to the free case, as Q passes 1/2, the most relevant angular momentum 
/ in the sum (13.701) changes from I = —1/2 to / = 1/2, and ai Q = 1/2 we have the OPE 
(I3.52P with two degenerate impurity operators. 

To find the nontrivial impurity exponent we need to solve eq. (13.671) for a{Q). We are, 
therefore, after the propagator D{x,x',Q) with x —>■ x'. We could, in principle proceed 
as in the free case. Namely, make our flux-tube uniform along D — 2 spatial dimensions 
(introducing a convergence factor v^~^ into (13.701) ). perform the integrals in (I3.70p with 
r = r' , T = 0, perform the sum over the angular momenta /, take 6^0 and D 3. 
However, unlike in the free case, the sums over angular momenta cannot be now performed 
analytically in terms of hypergeometric functions (with nice analytic continuation for ^ 0). 
The sum over / can still be performed numerically, however, the convergence is rather slow. 
Nevertheless, we have been able to determine a{Q) numerically using this method. However, 
this method is less suitable for finding the electric field coefficient C{Q), which requires us 
to differentiate the propagator at 6' = 0, making the convergence properties of the series 
even worse. 

Instead, we shall use a different method, generalizing the integral form of the propagator 



(I3.58P derived in Ref. |47| to the present problem. As shown in Appendix [Bl the twisted 



propagator at M = oo is given by. 



D{r,e) = — dutanh{nu)^==U^^{e) (3.73) 



with Vy{Q) still given by eq. fl339|) . 

Now, a{ff) can be determined from (I3.67p . 



= lim(D(r, Q) - D{r, 9,Q = 0)) (3.74) 

9 — »0 



V 



sinh(27rv^^ + a) sinh(27rz/) 



47rr 7o \\/v'^ + a cosh(27rVz/2 + a) - cos(27r(5) cosh(27rz/) - 1 

(3.75) 

Eq. (I3.74P can be solved numerically for a[Q\ However, before we do this, let's verify our 
claim that (iA) = (i.e. a = 0) is not sufficient to satisfy (I3.66P for finite Q. Indeed, from 
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-a(Q) Coefficient a(Q) of the Lagrange Multiplier 

0.15 - 
0.125 - 

0.1 : 

0.075 - 
0.05 

0.025 ; 

- - Q 

0.1 0.2 0.3 0.4 0.5 

FIG. 4: Coefficient a{Q) of the Lagrange multiplier (iA(x)), see eq. (|3.64p . in the M = oo gener- 
alization of the dual theory. 



0)) = ^{{V^V{x))q-{V^V{x))q=o) 
--^(l-2|g|)tan(7r|g|), |Q| < 1 

(3.76) 

where expectation values in the first line of (13.761) are computed in the free theory. The 
precise value of expression (I3.76P is not very important for our purposes (although it is 
curious to note that like many quantities in the free theory it is non-analytic in Q at Q = 0). 
What is important for us is that expression ( ]3.76p is negative. This means that the twisted 
boundary condition effectively creates a repulsive barrier, leading to a decrease in V'^V 
compared to untwisted theory. To compensate for this decrease in the strongly interacting 
theory, we need {iX{x)) to provide an attractive potential for V particles. Hence, we conclude 
that a{Q) < for Q finite. One may be concerned that the square roots in expressions (13.731) . 
(I3.74P are ambigious for a < and z/^ < \a\. However, it turns out that these expressions 
do not depend on our choice of the sign for the square root as long as it is consistent. The 
numerical solution for a{Q) is shown in Fig. |H We note that this solution agrees with the 
one obtained using the spectral form of propagator (I3.70p . 

One can also attempt to use eq. (I3.74p to find a series solution for a{Q) near Q = 0. It 



(I3.74P we obtain 

lim (-D(x, x', Q,a = 0) — D{x, x' , Q = 

47rr Jq cosh(27rz/) — cos{2ij:Q) 
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Delta Impurity Anomalous Dimension of V(x) 
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fT' . . . , . . . . , . . . . , . . . . ^ Q 

0.1 0.2 0.3 0.4 0.5 

FIG. 5: Impurity anomalous dimension A-^p of the monopole operator V{x), see Eq. p.48p . 
computed in the M = oo generalization of the dual theory. 



is easy to convince oneself that, 

a(g) ^ -Q\ Q^O (3.77) 

Unfortunately, the integrand in eq. f l3.74p is quite singular at z/ — > for a ^ 0, Q — 
0, so that a systematic series expansion beyond the leading order is not straight-forward. 
Nevertheless, we believe that such an expansion exists and a{Q) is an analytic function 
of Q near Q = 0. Assuming such analyticity and using charge conjugation symmetry, 
a{Q) = a{—Q), one obtains, a{Q) ~ —Q^ + c^Q^ for Q — >• 0. Here C4 is a positive constant 
as the integral (13.741) diverges for a < —Q^. 

Having found a{Q) we immediately obtain the impurity anomalous dimension of the 
operator V (given by eq. ( 13.72^ ). see Fig. [5l This anomalous dimension is no-longer the 
trivial value A^^p = \Q\ of the free theory (I3.49p . Given the leading behaviour of a{Q) as 
Q — > (I3.77P and assuming analyticity of a{Q) we conclude that ^Ymp ^i^^ ^^^o be analytic 
at Q = (as opposed to the situation in the free theory). Moreover, 

ALp~v^Q'> Q^O (3.78) 

Finally, we can now compute the coefficient of the electric field C{Q). For the M-field 
generalization of the dual theory, we define the magnetic field by the same equation (13.360 as 
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c(Q)/M Coefficient C(Q) of the Electric Field 

0.08 



, , 



0.06 



0.04 



0.02 




FIG. 6: Coefficient C{Q) of the electric field, see eq. (13.20j) . The doted and solid curves corresponds 
to the strongly interacting theory at M = oo and the free theory respectively. 



for M = 1 theory, that is we consider the current associated with the global U{1) symmetry. 



{-zF,{x)) = {-2mV^d,V^ix)) = -2m \imid;-dy){V^ix)V^iy)) = -2mM \imid;-dy)Dix,y) 

x—>y ^ ^ 

(3.79) 

Due to our normalization of the U{1) current, the electric field induced will be of order M. 
Now, differentiating D{r, 6) in (13.731) and taking the symmetric limit as ^ 0, 

o . ^ 1 P , w ^ sin(27rg) 

-ideD(r,0) = / c/i/z/tanh vri/ , ^ ^' 3.80 

4:nr Jo cosh(27rv/^;2:f^) _ cos(27rg) ^ 

Using the values of a{Q) found earlier (Fig. Hj) and evaluating the integral (I3.80p numerically 
we obtain the coefficient C{Q), shown in Fig. [H] (dotted curve). Fig. [Hlalso shows the value 
of C{Q) in the free theory (I3.56P for comparison (solid line). 

Alternatively, we can use (13.801) to expand C{Q) in a series in Q. Using the leading 
behaviour (13.771) . we find, 

c(g) ^ M + o(g3)^ , g^o (3.81) 

We see that the leading term in (13.811) agrees with the one, which would be obtained by 
perturbation theory in Q in the large M limit (I3.26p . (I3.34p . It is also interesting to compare 
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eq. fl3.8ip to assymptotic behaviour of C{Q) in the free theory (13.570 . We see that the 
leading term C{Q)/M ^ ttQ/8 in both cases is the same, however, the subleading terms are 
different. The first subleading term in the free theory is non-analytic ~ \Q\Q, as opposed 
to the strongly interacting theory's analytic 0{Q^). Thus, we have been able to verify that 
the leading non-analyticity of C{Q) in the free theory disappears in the interacting theory. 
We actually expect that the interacting theory cures itself of non-analyticities in Q at all 
orders in Q. 

Finally, let us discuss impurity anomalous dimensions of higher flux operators for 
rational Q = p/q, as these are actual physical observables on the QED side of the duality. 
Once we go from M = 1 dual theory to its large M counterpart, there are many possible 
generalizations of the F^(a;) operator. Indeed, we can form different SU{M) multiplets 
out of q instances of SU{M) fundamental Va{x). We expect that these multiplets will 
have different (impurity) anomalous dimensions for M finite. However, for M = cx) all of 
these operators will have degenerate (impurity) anomalous dimensions. We can consider, 
for instance, the completely symmetric representation Vg{x) = {Va{x))'^, where a is some 
fixed index (no summation over a). Then, for M = oo, 

{V^{x){V^{y))^) = q\{D{x-y)y (3.82) 

Hence, just as in the free case, the operator Vg{x) has the impurity OPE (I3.53P with the 
corresponding impurity anomalous dimension, 

ALp(g) = gA^p (3.83) 

IV. EASY PLANE THEORY FOR GENERAL N 

We now turn to the general case of the model S + iSimp with a global U(l)^~^ symmetry. 
The results of the previous section with = 1 can be rapidly generalized, and will lead to 
a quantitative result for the scaling dimension of the monopole/VBS operator V near the 
impurity. 

A. Duality in the Easy Plane Theory 

In this section, we consider a theory with flavours of spinon fields Za {N does not 
necessarily have to be large), 

L = + \{d, - iA,)zJ^ + f/(z„) (4.1) 

Here, \J is some potential with the global f/(l)^ symmetry under independent phase rota- 
tions of the Za fields. The singlet component of this symmetry is actually gauged by the 
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field A^, 

U{1) : ^ e^'^^h^, A^^A, + 8^9 (4.2) 
while the non-singlet components are true global symmetries of the theory, 

U{ir-' : - e^^'-'^-z^ (4.3) 

where t", a = 1..N — 1 are the generators of the U{1)^^^ symmetry satisfying, Yla^a ~ 0- 
We require U to have a symmetry under the permutation of labels of z^ fields. We choose U 
in such a fashion that in the "condensed" phase of the theory, it favours non-zero expectation 
values of all components of the Za field, so that the vacuum manifold of the theory is a torus, 
(S*^)^ (here we temporarily forget that the singlet symmetry is gauged). For N = 2 the 
theory under consideration is believed to describe the phase transition in the easy-plane 
ant if erromagnet . 



22 



We would like to dualize the theory (14.11) . Similar theories were dualized in Ref. 
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50, and here we will present a related discussion. An exact duality on the lattice 
appears in the compainon paper—, but we can write down the form of the dual action from 
very general considerations. Let us first identify the dual degrees of freedom. We go to the 
condensed phase of the theory (14.11) . where all (za) 7^ 0. Then, we can have vortices in any 
component of the z^ field. Formally, the homotopy group, 7ri((5'^)^) = Z^. So, we have N 
types of vortices, which become the degrees of freedom of the dual theory V^, a = 1..N. 

These vortices are global, rather than local. Indeed, let's consider a vortex in the first 
component Zi, 

zi{x) ^ ve'^^'^\ Zq, ~ v, a 7^ 1, |x| ^ 00 (4.4) 

where A(x) winds from to 27r as one goes around a contour out at infinity surrounding 
the vortex. Then, this vortex corresponds to a space-time dependent transformation of the 
vacuum (g^]), (O, with, 9{x) = ^A(f) and r(f)r = (1 - l/N,-l/N,... - l/N)X{x). 
Thus, our vortex possesses a winding both in the local and in the global symmetry group. 
The winding in the local U{1) group will be canceled by the gauge field, 

A,ix) = d.eix) = j^d.Xix) (4.5) 

hence our global vortices carry a magnetic fiux $ = 27t/N.— Therefore, under the flux 
symmetry (13.41) . the fields Va should transform as, 

VM ^ e2™(")/^V;(a;) (4.6) 

This fact will be crucial for the analysis to follow. 

The winding in the global group will lead to a long-range Coulombic interaction between 
our vortices. We will need dynamical gauge fields in the dual theory to give rise to this 
interaction. However, if we have a unit winding in each component of the z field, our vortex 
becomes completely local, and carries total fiux 2tt. We can think of such a local vortex as 
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a composite of global vortices of different types. The creation operator for this flux-tube, 
therefore, will be, 

V{x)=l[V^{x) (4.7) 

a 

Since the local vortex carries flux 27r, we can also associate the operator (14. 7p with the 
monopole operator of the direct theory. Indeed, given ( 14. 6p . under the flux symmetry ( 13. 4p . 

V{x) e2™Wv(x) (4.8) 

which is the correct transformation law for the monopole operator ( 13. 6p . 

We expect local vortices to interact by short range forces. Therefore, the operator (14. 7p 
should not be charged under the emergent gauge fields of the dual theory. 

We are now ready to write down the dual theory, 

i 

Here -B^ = , a = 1..N—1, are emergent dual gauge fields, which couple to the non-singlet 
currents. = e^y\dyB'^ are the corresponding field strengths. The dual potential UiVa) 
is chosen to have the same properties as the direct potential U: it has a U{1)^ symmetry 
under independent phase rotations of the fields Va and a symmetry under permutation of 
labels of Va fields. Moreover, it favours (V^) ^ for all a in the condensed phase of the 
dual theory. Thus, the theory (14.90 has a local U{\)^^^ symmetry, 

f/(l)^-i: \4(x)^e*^''(^)*^^\4(x), B;^B; + d,r (4-10) 

as well as the global U{1) flux symmetry of the direct theory ( 14. 6 p (which we have promoted 
to a local symmetry by introducing a non-dynamical source field H^. As required, the 
monopole operator (14.70 is invariant under the local U{1)^^^ symmetry of the dual theory 

dnoD. 

The theory (14.90 also has a global U{1)^~^ symmetry associated with conservation of 
fluxes of the A^ — 1 emergent gauge fields. This topological symmetry can be identified with 
the Noether's symmetry ( 14.30 of the direct theory. 

B. Wilson Loops in the Easy Plane Theory 

Now, we would like to apply the duality discussed in the previous sections to study the 
properties of Wilson loops in the U{1)^^^ symmetric theory (14. ip . Recall, that to represent 
Wilson loops we must use a source field given by (I3.1O0 . As discussed for the case of 
A^ = 1 theory, the effect of such a source field on the dual action ( 14. 9 p is to introduce a 
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twisted boundary condition for the vortex fields, 



V^{e = 27r) = e-'^^'^/JVy^^^ ^ g) (4.11) 

where Q is the charge of our Wilson line. The physical origin of the factor 1/A^ is the 
fractional charge 2tt/N of the vortex fields Va under the flux symmetry. Thus, we come 
to the amazing conclusion that the universal physics in the planar model is periodic in the 
charge Q of the Wilson line, with period Q = N. This is a generalization of the Q = 1 
periodicity of single flavour QED discussed before. As explained in a companion paper—, 
we expect that this Q N periodicity is a feature of the easy plane theory and does not 
generalize to the case with the full SU {N) invariance. 

Now, we would like to discuss more quantitative features of Wilson loops in the pla- 
nar model. In particular, we would like to find the impurity anomalous dimension of the 
monopole operator (13.481) and the coefficient of the electric field (I3.20p at the critical point 
of the theory. We note that as in the = 1 case, we can easily dualize the magnetic field 
by differentiating the dual action with respect to the source field H^, 

HF,) = ^-^{V^'dM (4.12) 

with D^Va = {d^ - iB'^ - ^H^)V^. 

To find Aj^p and C{Q), we follow the procedure established for the A^ = 1 case in section 
IIII Dl and perform a large M expansion of the dual theory (14. 9 p . Namely, we promote each 
field Va to an SU{M) multiplet, V^, i = 1..M. Moreover, we replace the soft potential UiVa) 
by a hard constraint, Yli iKvP — ^1 9i each a = 1..N. This constraint will be enforced 
by a set of A^ Lagrange multipliers Aq,. Thus, our Lagrangian becomes, 

L = J2 \{d, - ^B-^ - ^H,)V:\' + J2^Xa{K\' - -) (4.13) 

In f l4.13p we have also dropped the kinetic term for the gauge fields, as near the critical point 
such operators will be irrelevant. In addition to the U{1),^ global fiux symmetry and the 
U^l"jN-i iQf^oi symmetry of the original M = 1 action, the theory (14.130 also has a SU{M)^ 
global symmetry under independent SU{M) rotations of the A^ M-tuplets V^. We note that 
the various SU{M) multiplets talk to each other only through the gauge fields B°. 

We would like to generalize the observables of the M = 1 theory to the large M case. 
The magnetic field (I4.12p is generalized trivially, 

HF,) = ^-^m^^.V^) (4.14) 
The monopole operator (14. 7p on the other hand, now carries indices under the SU{M)'^ 
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group, 

V{x),,...^ = l[v:-{x) (4.15) 

a 

The insertion of the Wilson loop source if^ is again equivalent to the twisted boundary 
condition (14. lip . 

We now perform a large M expansion of the theory (14.131) with the twisted boundary 
condition (14.111) . keeping fixed. We will be only able to make computations for M = oo. 
We are interested in the physics at the critical point. We expand the theory about the 
saddle point = (this is a saddle point as the twisted boundary condition (14.111) does 
not couple to the non-singlet sectors of the theory^). As usual, the fluctuations of these 
gauge fields about the saddle point will be suppressed by powers of 1/M. Thus, at M = oo, 
we are left with decoupled instances of the Lagrangian (I3.60p that has been discussed at 
length for the case of A^ = 1 theory. The only difference is the replacement, Q ^ Q/N in 
the boundary condition (13.191) . Hence, we conclude, 

{V{x),,..,VHx\.,^) *=°° n(^"(^)(^i")H^')) = D{x,y,Q/Nfl[6,^,^ (4.16) 

a a 

where D{x,x' ,Q) is the propagator in the A^ = 1 theory (I3.60p with the twisted boundary 
condition (13.191) at M = oo. The asymptotic behaviour of this propagator for r -C r' is 
given in eq. (I3.7ip . Thus, the asymptotic behaviour of the correlation function (I4.16P for 
r -C r' is 

(V(x),...,Vt(x'),,.,,) ^ ^^yViQ/N^+^iQ/N) e-O^G(r/r') H'^w., \Q/N\ < 1/2 

(4.17) 

where G is some (known) function. Hence, the monopole operator V{x) in the planar A^ 
component theory has the impurity OPE, 

V(x, r) ~ Ixl^^^-pe"^^^ Vimp(r) for \x\ (4.18) 

with 

A^^p = N^iQ/Ny + a{Q/N) = NAl,^,{Q/N), \Q/N\ < 1/2 (4.19) 

where the monopole impurity anomalous dimenension A^^-^^{Q) in the A^ = 1, M = cxo 
theory is given by Fig. [5l 

From OPE (I4.18p . we observe that for integer Q the monopole operator is single valued 
under 9 9 + 2n, even though the dynamical fields of the theory obey twisted boundary 
conditions (14. lip . We also note that formulas (I4.18P and (I4.19P are correct only for \Q/N\ < 
1/2; for other values of Q they should be extended by periodicity Q Q + N. 

We can now take the A^ oo, Q-fixed limit of (I4.19p . Using the assymptotic behaviour 
(I3.78p . Af^p ~ Q'^/N. Thus, the impurity anomalous dimension of the monopole operator 
is of order 0{1/N) for A^ — * oo in the easy plane theory. It is interesting to note that, as 
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discussed in section |TT1 this is also true of the theory with a full SU{N) symmetry. At this 
point, it is not clear whether this is just a coincidence. 

Finally, let us discuss the special point Q/N = 1/2. Our interest in this point is not 
purely academic, as we expect iV = 2, Q = 1 to correspond to the physical case of a single 
impurity in an easy plane antiferromagnet or superfluid. We recall that at this point the 
propagator D{r, r', 9, r) for r -C r' is dominated by two angular momenta, / = ±1/2, 

+ (4.20) 



So that 



« ( i^'^'f: (4.21) 



m=0 

Hence, using (14.161) . the correlation function of two monopole operators is dominated by 
angular momenta / = —Q, —Q + 1..Q — 1,Q for r <^ r'. So, we conjecture the operator 
product expansion, 

Q 

V(f,r) ~ J2 ci\x\^^ e~'^^ Vi{t) for \x\ (4.22) 

l=-Q 

At M = oo all the operators Vi have degenerate impurity anomalous dimensions . As 
discussed in section IIII CI the anomalous dimensions of operators with oposite angular mo- 
menta are equal by CP symmetry emergent at the Q/N = 1/2 point. However, there is no 
fundamental reason why anomalous dimensions of operators with different values of / should 
be equal. Thus, we expect the degeneracy to be lifted at higher orders in 1/M expansion. 
Therefore, unfortunately, the question of whether the OPE (14.221) will be dominated by / = 
or by finite / is beyond the reach of our calculation. Nevertheless, our calculation ai M = oo 
predicts for the physically relevant case of A^ = 2, Q = 1, 

N = 2,Q = 1 (4.23) 

The emergent CP symmetry at the point Q/N = 1/2 means that quantum fiuctuations 
manage to render the states of Figs. [2] and Fig. [3] degenerate in the long-wavelength limit. 
We remind the reader the CP symmetry is due to the emergent Q ^ N periodicity of the 
easy plane theory. No such periodicity is expected to occur in the full SU{N) symmetric 
theory, where the impurity OPE is dominated by a single operator with a definite angular 
momentum as in eq. (11.121) . 

For completeness sake, we also discuss the coefficient C{Q) of the electric field. From eq. 
(I4.14P at M = oo we obtain, 

C{Q)=Cn=i{Q/N) (4.24) 
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where the coefficient Cn=i{Q) in the = 1, M = oo theory is given by Fig. [61 We note 
that for Q/N = 1/2 the electric field vanishes, as it should, by the emergent CP symmetry. 

V. CONCLUSIONS 

This paper began with the theory S in Eq. (II. ip for square lattice quantum antiferro- 
magnets in the vicinity of a Neel-VBS quantum phase transitions. We considered generic 
local deformations of the antiferromagnet, and argued that they could be classified into two 
categories. The first category, illustrated in Fig. [1], is a modulated exchange impurity: we 
found an enhancement of VBS order, characterized by the exponent in Eq. II. 9[ The second 
category was realized by a missing or additional spin {e.g. Zn or Ni impurities on Cu sites), 
shown in Fig. [2l For this case we found that VBS order was suppressed by the appearance 
of a VBS pinwheel, as in Fig. [2l and characterized by the scaling properties discussed in 
Section II Bl 

The results of this paper should be useful in numerical studies of the quantum phase 
transition between the Neel and VBS state^*^. By enhancing an exchange constant as in 
Fig. [H and measuring the decay of the average VBS order parameter away from the impurity, 
the exponent can be estimated from Eqs. (11. 6111. 9p . There will be no mean VBS order 
in the vicinity of a missing spin impurity as in Fig. O However, the spatial dependence in 
the VBS susceptibility is fixed by Aj^p in Eq. fll.l2p . The positive value of A?^p indicates 
that the VBS susceptibility should be suppressed near such an impurity. 

In STM studies of the cuprates, we have noted earlier the demonstration of bond-centered 
charge order in the local density of states by Kohsaka et al.—. A numerical analysis of the 
pinning of such charge order by modulated exchange impurities (in the class in Section II Ap 
has also been carried out^"^. However, it is also experimentally possible to induce "missing 
spin" impurities (in the class of Section II BP by replacing the Cu sites with Zn and Ni impu- 
rities. There have been STM studies of such impurities^i^i^, and it be of great interest to 
carefully examine the nature of the bond-centered modulations in the vicinity of such impu- 
rities. If we assume that the "stripe" instability is primarily associated with the appearance 
of magnetic order— then the theory of the enhancement of magnetic order near such 
impurities^'^i should apply: we should therefore expect an increase in the strength of the 
density of states modulations in this model. In contrast, if we assume a VBS theory of 
the modulations, then in the impurity model of Section II B| the bond-centered modulations 
should be suppressed. The experimental situation could well include both effects, compli- 
cating the interpretation. However, evidence for VBS pinwheel configurations like those in 
Fig. [2] would lend strong support to the VBS theory. 
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APPENDIX A: ELECTRIC FIELD IN THE FREE THEORY 



We wish to use eq. fl3.55p to compute the electric field in the free theory. We start 
from the equation for the propagator f l3.38p and promote the uj integral to run over D — 2 
dimensions, as discussed in section IIII C[ obtaining 

^(^'^) = ^Ip^EC^ / duu--^,iurr (Al) 

^ r((3-D)/2) 1 .^e-'^ r(|/|+D/2-l) 

(47r)(^-i)/2 27r r(|/| - D/2 + 2) ^ ^ 

We see that the pref actor diverges for D = 3. However, at D = 3 the sum over angular 
momenta becomes ^ = S{6). So, we have to first perform the sum over angular momenta 
and then take the limit D ^ 3. The sum over angular momenta can be performed in terms 
of hypergeometric functions, giving for < Q < 1, 

p^,m _ r((3-/^)/2) e-QV , r(D/2-Q) , 



(A3) 



where F denotes the Barnes extended hypergeometric function. One can check that for 
D = 3 the expression in brackets in flA3l) vanishes, cancelling the pole in the prefactor. 
Now, differentiating with respect to 6*, 

■f^jii (^\ r((3-D)/2) e--QV (i-g)r(D/2-g) 

-^deD{r,e) = [ r(3-D/2-g) " Hih D/2 - Q}, {3 - D/2 - Q} 

+ r(f-^D/2lg) ^"'^^^^' /^/2 - g + 1}, {4 - D/2 - Q}, e^^) 
QTiD/2 - 1 + g) ^^^^^ Z)/2 - 1 + g}, {2 - D/2 + Q}, e^^) 



r(2-z}/2 + g) 
r(D/2 + g) 



r(3-D/2 + g) 



e-^F({2, D/2 + g}, {3 - D/2 + Q}, e^'')) 



According to (13.551) . to compute the electric field we need to take the limit as ^ of 
(1A4I) . Strictly speaking this limit does not exist as the hypergeometric functions blow up as 
6^0 (that is when the last argument goes to 1). However, we note that only the imaginary 
part of (1A4I) becomes infinite as 6' — > 0, while the real part has a well-defined limit. The 
expectation value of electric field {—iEr) = —{—iFg) should be real. Thus, we can drop 
the infinite imaginary part. Moreover, the imaginary part is antisymetric under 9 — > —9, 
so the "symmetrized" limit of (IA4p exists. It turns out that this symmetrized limit can be 
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obtained by the formal summation formulas, 

F({l,a},{6},l) = (A5) 
a — + 1 

So, taking 0, plugging flA5l) into flA4l) and performing a few manipulations, 

_ -^r,{f) = n ^_ (2Q-l)r((l-/^)/2)r(D/2 + g-l) / sin(7r(D/2 + g)) \ 1 
' ^ u,rj 2^+27r(^+i)/2 r(l-D/2 + g) Vsin(7r(Z}/2-Q)) 

(A7) 

Taking the limit D ^ 3, 

- tdeD{e = 0,r) = --^{2Q - l)Han{7rQ) (A8) 

o/vrr 

Finally, plugging into (13.551) we recover (13.201) with 

C{Q) = l{l-2Qftan{nQ), 0<Q<1 (A9) 

o 

We remind the reader that all the manipulations above have been performed for < Q < 1- 
The function C{Q) can then be extended to other values of Q by periodicity. In particular, 
extending to the range |Q| < 1, 

C{Q) = 1(1- 2\Q\YtaninQ), \Q\ < 1 (AlO) 



APPENDIX B: INTEGRAL FORM OF THE TWISTED PROPAGATOR 

In this section we review the derivation of the integral form of the twisted propagator 



(I3.58P given in Ref. |47l. We use this integral form to compute the electric field (13.551) and 
show that it is in agreement with the result obtained using spectral representation of the 
propagator (see Appendix [A|) . We also indicate how the free twisted propagator should be 
modified in the strongly interacting M = oo theory. 

Recall the free massive propagator in 2D (without any twisted b.c.) obeys, 

i-d'^ +m'^)D{x,x) = 6{x-x') (Bl) 

and is given by, 

1 1 /"°° 

D2{x,x') = —Ko{m\x-x'\) = —r duKUmr)KUmr')e'"'e-''\'-''\ (B2) 
2vr 2tt^ 

where the integral representation is valid for \6 — 6'\ < 27i. The BesselK functions of 
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imaginary argument satisfy the equation, 



-li (^i) - ^ + ) KU^^) = (B3) 



r dr \ dr J 

Hence the functions Kiy{rnr)e^^^ are in the kernel of the operator — 9| +m 
'^'M^ + Applying this operator to D2{x, x') we learn, 



r or ^ or ' 



I poo 

duuKi^{mr)Ki^{mr')e'"' = -6{r - r') (B4) 
J-oo r 

This identity will be useful to us later. 

Now, we want to modify the propagator ( ]B2p in such a way that it satisfies the twisted 
boundary conditions ( 13.190 . Let's first symmetrize equation (1B2I) with respect to u by noting 
Kii, = K_i^. Then, 

1 /""^ 

D2(r, r', d-d') = — diyKUmr)KUmr') cosh(i/(7r - - ^D) (B5) 
Now, we can generalize, 

1 r°° 

L>2(r, r', e,Q) = — dpK,,{mr)K,,{mr') sinh(7rz/)f/,(0) (B6) 

^TT^ J-oo 

where 

smh(7rz/j 

D2{r,r' ,9,Q) still satisfies eq. fIBip since, as noted above, the functions Ki^{mr)e'^'^^ are 
in the kernel of — 9| + m^. It remains to find c(i/) such that the propagator (1B6P obeys 
boundary conditions fl3.19p . After a few manipulations one arrives at, 

e-'^''iQ^s^(o) smh{u\9\) + sinh(z/(27r - \9\)) 
^ ~ cosh(27rz/) - cos(27rg) ^ 

Next, one uses the identity, 

TT /""^ 1 1 

^h.[T:v)Kiy[mr)Kij,{mr') = — i rfw Jo("^(2rr') 2 (cosh(M) — cosh^2)^) sin(z/u) (B9) 



sml 



where ^2 > is defined by. 



2 I /2 

^ _i_ .y^' 

co«h6 = ^^ (BIO) 



Substituting this into flB6p , 

2 roo POO 

D2ir,r\e,Q) = — du Jo(m{2rr')^ {cosh{u) - cosh ^2)^) du U^{e) sm{uu) (Bll) 
if, Jo 
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We are mostly interested in the propagator with r = r\ 

POO roo 

D2{r = r',e,Q) = — du Jo(mrV2{coshu - l)^) duU^{e) sm{uu) (B12) 
2vr Jo Jo 



In principle, it is possible to perform the integral over i/ analytically in (IBlip (see Ref. 1471 ) . 
This, however, will not be very benificial for our purposes. Instead, let's proceed directly to 
the three-dimensional massless propagator, obtained by integrating over the mass parameter 
of the two dimensional propagator (13.381) , 

Dir,9) = -—^ du—- / duUM^H^u) (B13) 

ivr^rv^ Jo (cosh -u — 1)2 Jo 

where we have computed only the 3 dimensional propagator with r = r', r = t'. Now, 
performing the integral over u, 



D{r,e) = / diytaBh{TTu)U^{e) (B14) 

4vrr Jq 

To find the electric field we again use eq. (I3.55p . 

1 P, sin(27rg) cosh(z/0) 



idgD(r,9) = —- — / (iz/ z/tanh(7rz/) 
Anr Jq V 



cosh(27rz/) — cos(27r(5) 



(cos(2.g)oosh(.g)-ooshH2.-|»l)))) 

^ ^ ^ cosh(27rz/) - cos(27rg) ^ ^ ^ 

Again, the real part of —ideD{r,9) has a well-defined limit as ^ — > 0, while the imaginary 
part is antisymmetric under 9 —9 and diverges as ^ — > 0. So the "symmetrized" limit is 
given by, 

^ , , 1 [°° , sin(27rQ) tanh(7ri/) 1 / , 

-m,D(r,B = 0) = -- 1 Jp^) = -^(21(31 - l)Han(.Q) 

(B16) 

in agreement with an earlier computation (\A8\i based on spectral decomposition. Thus, 
C{Q) is again given by expression (13.561) . 

Now we generalize the above derivation of the twisted propagator to the strongly inter- 
acting M = oo theory. The strongly interacting theory differs from the free theory by the 
aditional space- varying potential (iA(x, r)), so that the propagator satisfies, 

i-d' + j^)D{x, x', Q) = 5{x - x') (B17) 

We again rewrite D{x,x' ,Q) in terms of the two dimensional massive propagator 
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D2{x, x', rn^, Q) as in eq. fl3.38p . The two dimensional propagator satisfies, 

We need to generalize the two-dimensional, massive, twisted, free propagator flB6P so that 
it obeys the above equation. We observe that the function Uy{6) (]B7p . (IBSP satisfies, 

1^ = u'UM - 2^5(^) (B19) 
Now combining eqs. (1B3|) . (]B4p and (IBlOp . we find that. 



D^ir, r', e,Q) = — duKUmr)KUmr') sinh(7rz/)^^=f/^^(^) (B20) 
Jo + a 

satisfies (IBISP as needed. Proceeding as above from two to three dimensional propagator, 
and setting r = r\ t = t' 

1 V 
^^"^^^^ = 1^ rfz/tanhM^==f/^(^) (B21) 
47rr Jq + a 
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